Solitonic fullerene structures in light atomic nuclei 
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The Skyrme model is a classical field theory which has topological soliton solutions. These solitons 
are candidates for describing nuclei, with an identification between the numbers of solitons and 
nucleons. We have computed numerically, using two different minimization algorithms, minimum 
energy configurations for up to 22 solitons. We find, remarkably, that the solutions for seven or more 
solitons have nucleon density isosurfaces in the form of polyhedra made of hexagons and pentagons. 
Precisely these structures arise, though at the much larger molecular scale, in the chemistry of 
carbon shells, where they are known as fullerenes. 
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PACS Numbers : 

The Skyrme model jjj was first proposed in the early 
sixties as a model for the strong interactions of hadrons, 
but it was set aside after the advent of quantum chro- 
modynamics (QCD). Much later Witten 0] showed that 
it could arise as an effective description at low energies 
in the limit where the number of quark colours is large. 
Subsequent work || demonstrated that the single soliton 
solution (known as a Skyrmion) reproduced the proper- 
ties of a nucleon to within an accuracy of around 30%; 
quite an achievement, given that there is, at present, no 
practical way of calculating the properties of nuclei from 
QCD via, for example, lattice gauge theory. 

In order to study nuclei of larger atomic number one 
first needs to compute the minimal energy configurations 
of multi-solitons, since in the Skyrme model there is an 
identification between the numbers of solitons and nucle- 
ons. Here, we present the results of an extensive set of 
simulations using two very different approaches designed 
to compute the minimal energy solutions for upto 22 soli- 
tons. With a small number of caveats, these results estab- 
lish an attractive analogy with fullerene cages familiar in 
carbon chemistry jJH . Although these classical solutions 
must first be quantized before a final comparison with 
experimental data can be performed it is expected that 
quantum corrections will be relatively small, since we are 
dealing with solitons, and so the classical solutions will 
contain important physically relevant information about 
the properties of nuclei. 

The Skyrme model is defined in terms of an SU(2) 
valued field E/(x), with an associated static energy 
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Note that the two physically relevant constants which 
would appear in front of each of the two terms in the 
most general version of (Q) have, for convenience, been 
scaled out by an appropriate choice of energy and length 
units. For finite energy, we impose the boundary condi- 
tion U(oo) = I, and pions are described by the usual 



quantum field theory treatment of fluctuations of the 
Skyrme field around this vacuum value, but nucleons 
arise in a very different manner, as classical soliton solu- 
tions. 

The boundary condition implies a compactification of 
the domain, and therefore U is a map from compacti- 
fied 1R 3 ~ S 3 h- > S 3 , since S 3 is the manifold of the 
target space, the group SU(2). Such mappings have non- 
trivial homotopy classes characterized by an integer val- 
ued winding number, which has the explicit representa- 
tion 
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B, which stands for baryon number or the number of nu- 
cleons, is often referred to as the topological charge and 
is the number of solitons in a given field configuration. A 
simple manipulation of equations (^) and (|2|) allows one 
to deduce the Faddeev-Bogomolny (FB) bound E > \B\. 
Generically, a charge B field will have an energy density 
E (the integrand of ([[])) and a baryon density B (the inte- 
grand of (H)) both of which consists of B well-separated 
lumps localized in space. However, as we discuss below, 
this is not the case for the minimal energy fields in which 
the solitons are close together. 

The mathematical problem is to find, for each integer 
B, the field U which minimizes the energy (Q) subject 
to the constraint (|^). This can be addressed by numeri- 
cal algorithms designed to minimize cither a discretized 
version of the energy or, equivalently, by solving a dis- 
cretized version of the second-order field equations which 
follow from the variation of (|l|) . This first approach is a 
very demanding computational exercise (see ref. || for a 
detailed discussion), requiring the use of modern parallel 
supercomputers, but results for low charge {B < 8) were 
found using this method J7]— p~0| . The results presented 
here extend this numerical approach up to B — 22. In ad- 
dition, we have applied a second, very different, technique 
to the construction of minimal energy solitons which not 
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only allows us to have greater confidence that the numer- 
ical solutions we have constructed are indeed the global 
minima, but in addition provides a good analytic approx- 
imation to the numerical solutions, making it much easier 
to identify their structure and symmetries. 

Our second approach makes use of the remarkable fact 
that minimal energy Skyrmions can be approximated 
by an ansatz involving rational maps between Riemann 
spheres p| ; a result which we will further confirm. The 
Skyrme field is a map U : R 3 i— > S* 3 , so it is not imme- 
diately obvious how to obtain such a map from rational 
maps which are between spheres S 2 i— > S 2 . Briefly, the 
domain S 2 of the rational map is identified with concen- 
tric spheres in ]R 3 , and the target S 2 with spheres of 
latitude on S 3 . To present the ansatz it is convenient to 
use spherical coordinates in ]R 3 , so that a point x S R 3 
is given by a pair (r, z), where r = |x| is the distance 
from the origin, and z is a Riemann sphere coordinate, 
namely z = tan[0/2] exp[i0] where 9 and <j> are the normal 
spherical polar coordinates. 

Now, let R{z) be a degree B rational map, that is, 
R = p/q where p and q are polynomials in z such that 
max[deg(p), deg(q)] — B, with no common factors. Given 
such a rational map the ansatz for the Skyrme field is 



U (r, z) = exp 



1 + \R\ 2 



-\R\ 2 
2R 



2R 

\R\ 2 - 



(3) 



where f(r) is a real profile function satisfying the bound- 
ary conditions /(0) = tt and /(oo) = 0. This is deter- 
mined by minimization of the Skyrme energy of the field 
(|J) given a particular rational map R. The ansatz yields 
an exact solution for B — 1 and it was shown in ref. 
jll| that for 2 < B < 8, suitable maps exist for which 
the field (^) is a good approximation to the numerically 
computed solutions, in the sense that the symmetry is 
identical and the energy is only one or two percent above 
the numerically computed values. 

Substituting the ansatz (||) into the energy (|l|) pro- 
duces an energy function on the space of rational maps, 
which we denote by I(R), given by 
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Therefore, our second approach to computing minimal 
energy Skyrmions is to search the (finite dimensional) 
parameter space of general degree B rational maps to 
find the one which minimizes 2"(i?) , using a powerful nu- 
merical minimization technique known as simulated an- 
nealing pm . 

Clearly, this procedure is not guaranteed to find the 
minimum energy Skyrmion since the topography of the 
rational map space may be slightly different to that of 
the full non-linear field theory, but as we shall see for 
the most part it works well, only encountering difficulties 
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FIG. 1. The baryon density isosurfaces for the solutions 
which we have identified as the minima for 7 < B < 22, and 
the associated polyhedral models. The isosurfaces correspond 
to B = 0.035 and are presented to scale, whereas the polyhe- 
dra are not to scale. 
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when there are two or more Skyrmion solutions, either 
saddle points or genuine local minima, with very simi- 
lar energies. We use the first minimization technique as 
a check, and in the small number of cases where relax- 
ing well-separated clusters consistently yields a different 
solution for a wide range of initial conditions, the sym- 
metry of the Skyrmion solution is identified by eye from 
the baryon density isosurface, and an approximate ratio- 
nal map can then be found by relaxing in the rational 
map space restricted to have the correct symmetry. In 
such cases the values of X for the different solutions are 
usually very close ||. 

The results of applying the two minimization tech- 
niques in this way are presented in table | for 1 < B < 22 
and pictorially in Fig. |l| for 7 < B < 22. In all but a 
small number of cases (B = 10, 14, 16, 22) we find that 
the minimum energy Skyrmion and that in the rational 
map space coincide. Furthermore, for each of these spe- 
cial cases, except B — 14, we were able to find a map 
with the same symmetry. For B = 14 we were prevented 
from finding a rational map approximation for the true 
minimum since its symmetry group, C2, is a subgroup of 
that of the minimum energy rational map, D^. 

The baryon density isosurface can be associated with 
a polyhedron whose edges and vertices coincide with the 
regions in which the baryon density is localized. Exami- 
nation of the solutions shows that, with the exception of 
B = 9 and 13, the associated polyhedra are trivalent with 
A(B — 2) vertices (the Geometric Energy Minimization 
(GEM) rules) as predicted in ref. 0, and for B > 7 they 
comprise of 12 pentagons and 2(B — 7) hexagons. Such 
structures are common in a wide range of physical appli- 
cations, and have become a hot topic in carbon chemistry 
where they correspond to shells with carbon atoms placed 
at the vertices, the most famous being the icosahedrally 
symmetric Buckminsterfullerene Cqo structure, which is 
also the traditional soccer ball design. For this reason we 
shall refer to such solutions as being of the fullcrcnc type, 
with the prediction, spectacularly confirmed by our re- 
sults in all cases except B — 9 and B — 13, that the poly- 
hedron associated with the Skyrmion of charge B, has a 
structure from the family of carbon cages for 

We had predicted in ref. (l(J that the Buckyball Cqo 
configuration would be found for B = 17 and indeed 
an approximate rational map description was found in 
ref. |fj"lf . Here, we see that such a solution is the mini- 
mum energy solution of the full non- linear field equations 
and in the rational map space. We see also that a large 
number of the other solutions have platonic symmetries 
which, from the mechanical point of view, implies the 
structure packs well. It would appear, therefore, that 
such structures may be preferred over less symmetric ones 
in the minimization procedure. We should note, however, 
that this is not always the case and rational maps with 
platonic symmetries can easily be found, for example at 
B = 9, which do not give minima 0. 
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TABLE I. A summary of the symmetries and energies of 
the Skyrmion configurations which we have identified as the 
minima. Included is the ionization energy (Ib) — that re- 
quired to remove one Skyrmion — and the binding energy 
per Skyrmion (AE/B) — that required to split the charge B 
Skyrmion into B charge one Skyrmions divided by the total 
number. 



The polyhedra found for B = 9 and 13 do not obey 
the GEM rules, nor are they of the fullerene type, since 
they contain four-valent links. They can, however, be 
related to a fullerene via the concept of symmetry en- 
hancement, as follows. A very common structure within 
the fullerene polyhedra is two pentagons separated by 
two hexagons. If the edge which is common to the two 
hexagons is shrunk to have zero length, the four poly- 
gons then form a C4 symmetric configuration containing 
a four valent bond. For the case of B = 9, the polyhedron 
can be thought of as being created from a D2 symmetric 
fullerene by the action of two such operations, and in the 
B = 13 case six operations can be used to convert another 
D2 configuration into one with O symmetry. Empirically, 
we see that each symmetry enhancement operation ap- 
pears to be accompanied by an equivalent one antipodally 
placed on the polyhedron, and single operations appear 
not to occur. 

We have computed the energies of the solutions which 
are presented in table [l| using the rational map ansatz 
to create initial conditions which are then relaxed un- 
der the action of the full non-linear field equations. It 
should be noted that these values are (for B > 1) al- 
ways a little less than the corresponding values computed 
solely within the rational map ansatz. On the discrete 
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FIG. 2. On top the ionization energy Ib plotted against B. 
Notice that the most stable solutions are those with the most 
symmetry, B = 4, 7, 17, while the least stable are those with 
little symmetry B = 5, 8, 14, 18. On the bottom the binding 
energy per baryon AE/B plotted against B. We see that for 
large B the binding energy appears to level out at around 
0.15 — 0.16 as one might expect in a simple model of nuclei. 

grid the computed value of the baryon number, £?dis, is 
less than the corresponding integer B suggesting that the 
finite difference approximations used to compute the en- 
ergy £?dis will under-estimate the true energy. Moreover, 
in the initial conditions one must impose the boundary 
condition U = 1 at the edge of the box. By using a 
wide range of different grid sizes and spacing we have 
shown H that the value of -Edis/i3dis can be computed 
accurately, and hence so can the true energy using the 
formula Eb = B * (E^i s / Bd; s )- We claim that our de- 
terminations of Edis/Bdis are accurate in the absolute 
sense to within ±0.001, and that the relative values are 
probably even more accurate. 

We have also computed the ionization energy Ib = 
Eb-i + Ei — Eb, which is the energy required to remove 
a single Skyrmion, and the binding energy per nucleon 
AE/B = Ei — (E/B) which is the energy required to 
separate the solution into B well-separated Skyrmions 
divided by the total baryon number. These values are 
tabulated in table Q and are plotted against B in Fig. |^. 
The ionization energy is largest for the most symmetrical 
solutions B — 4, 7 and 17, and is least for those with little 
symmetry, B — 5, 8, 14 and 18, which is very much as one 
would expect. The binding energy appears to increase to 
an asymptotic value of around 0.15 — 0.16. This is a clear 
consequence of the FB bound since it is linearly related 
to E/B. 



In fact the value of E/B appears to have an asymp- 
totic value which is around 6% — 7% above the FB bound, 
compatible with the value obtained for a hexagonal lat- 
tice [|l3| which is the limit of an infinitely large fullerene 
(the analogue of graphite in carbon chemistry). It is 
clear that an infinitely large shell is physically unlikely 
and that there probably exists a value B* such that for 
B > B* the solutions no longer look like fullerene shells. 
In such a case the solutions are likely to begin to look 
more like portions cut from the infinite Skyrme crys- 
tal Q whose E/B is only 4% above the FB bound. 
Another possibility is an intermediate state comprising 
of multiple shells [15| , although all the known configura- 
tions of this kind have much larger values of E/B. We 
have definitely shown that £?* > 22 and we believe that 
the connection between Skyrmions, fullcrcncs and ratio- 
nal maps will continue for much larger values of B. 
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